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Abstract. Models of galactic chemical evolution (CEMs) show that the shape of the stellar 
initial mass function (IMF) and other assumptions regarding star formation affect the resultant 
abundance gradients in models of late-type galaxies. Furthermore, intermediate mass (IM) stars 
undeniably play an important role in the buildup of nitrogen abundances in galaxies. Here I 
specifically discuss the nitrogen contribution from IM/AGB stars and how it affects the N/O- 
gradient. For this purpose I have modelled the chemical evolution of a few nearby disc galaxies 
using different IMFs and star formation prescriptions. It is demonstrated that N/O-gradients 
may be used to constrain the nitrogen contribution from IM/AGB-stars. 
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1. Introduction 

Much progress has been made in modelling nucleosynthesis in stars in general, but the 
origin of nitrogen still remains somewhat uncertain. In the classical picture, the apparent 
increase in nitrogen production with metallicity, as inferred from spectroscopy of galactic 
and extragalactic HH-regions, is a manifestation of nitrogen being either primary or 
secondary (produced from primordial elements or with heavier elements as "seeds"). 
However, reality often turns out to be more complex than our simplest models of it. The 
origin of nitrogen is probably no exception. 

During the last few years, reasonably accurate stellar nitrogen abundances have become 
available also for Population II stars (see, e.g., |Cayrel et al. 2004 lEcuvillon et al. 20041 



llsraelian et al. 20 04, Spite et al. 2005[ ). The rather large nitrogen abundances relative to 



iron seen in halo stars in the Milky Way suggest that a relatively large fraction of the 
nitrogen should come from Type II supernovae, especially at low metallicity. However, 
recent stellar evolution models, including nucleosynthesis, all appear to predict significant 
nitrogen production in intermediate-mass (IM) stars ([Marigo 2001, Izza rd et al. 2 004. 
Gavilan et al. 2005, Karakas & Lattanzio 2004). The typical production factor would be 
of the order ~ 10~ 3 of the initial stellar mass. If IM-stars of all metallicities produce that 
much nitrogen, it would require that metal-poor galaxies (or regions within a galaxy) 
with low N/O-ratios have very young stellar populations, in fact, they have to be so 
young that IM-stars have not yet made any significant contribution to the interstellar 
nitrogen abundance of these objects, unless the production of oxygen is much greater at 
low metallicity. 

In this work I have considered radial gradients of O/H and N/O in the Milky Way 
and five other nearby late-type galaxies in order to find additional clues to the origin of 
nitrogen. The advantages of considering the abundance gradients within galaxies, rather 
than over-all trends obtained from the global properties of galaxies, are that the same 
stellar nucleosynthesis prescriptions (yields) must be able to simultaneously reproduce 
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Table 1. Model parameters and adopted distances. See Sect. [3] for further details. 
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the radial abundance gradients in galaxies with quite different properties (sizes, total 
masses, gas distributions, etc.) and the very same yields must also be consistent with the 
fact that galaxies with similar O/H-gradients may have quite different N/O-gradients. 
Models of chemical evolution (CEMs) may thus be able to provide insight and constraints 
on stellar yields when several galaxies are modelled and analysed in parallel. 



2. Data 

The oxygen and nitrogen abundances are taken from Pilyugin et al. (2003, 2004). 
which provide a large compilation of data that is consistent in the sense that all abun- 
dances are derived using the same strong-line calibration (P-method, see |Pilyugin 2001a] 
|Pilyugin 2001bp . Data for HI and H2 (CO) were collected from the literature (see Table 
Q] for references) . To obtain a consistent and, possibly, a more correct set of data for the 
gas distributions, the conversion from CO-flux to H2 surface density has been redone 
using the oxygen-dependent calibration by Wilson (1995). 



3. Models 

I present detailed, numerical CEMs for the Milky way and five other nearby spirals, 
computed using the code described in detail in Mattsson (2008). Since the halo con- 
tributes very little to the present-day abundance gradients, except for at very large 
galactic radii, where the disc density becomes comparable with the halo density (see, 
e.g.,[Uhiappini et al. 2001] ), the model is constructed as a single exponential disc, formed 
by exponentially decaying baryonic infall, i.e., 



^inf 



M) 



r(r) 



exp 



to 



exp 



r(r) 



(3.1) 



where r is the infall-time scale and to denotes present time. The star formation rate is 
prescribed by a modified "Silk-law" ( |Wyse fc Silk 1989] IBossier et al. 20031) . i-e-, 



£*(r, t) = x] £* iM w(-R©, £0) ■ 



O (r) 
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(3.2) 



^0,Mw(-R©) 

where fio is the present-day (t = to) angular frequency of the disc (taken from Sof ue et al. 1 999) 
and 77 is a parameter regulating the star-formation efficiency relative to the Milky Way. 
To simulate the inside-out formation of a galactic disc, the infall time scale is assumed to 
have a linear radial dependence, r(r) = max(0, To + n r/r^), where to,ti are arbitrary 
constants and is the disc scale length. For tq = —1.3, T\ = 2.9 and = 2.8 kpc, the 
above relation is similar to that used by Romano et al. (2000) and Boissier & Prantzos 
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Figure 1. Radial trends for O/H (upper panels), N/O (lower panels) and the surface density 
of hydrogen. Black dots show the abundances from Pilyugin et al. (2003, 2004), red lines show 
models with yields from Gavilan et al. (2005, 2007) and blue lines show models with yields from 
van den Hoek & Groenewegen (1997). The solar symbols show the solar values according to 
Asplund et al. (2005). 

(1999). For the gas-rich galaxies (NGC 598, NGC 2403, NGC 6946), the infall rate is 
assumed to be constant at all radii, indicated by r — > oo in Table [TJ 

The nucleosynthesis prescriptions were taken from Chieffi & Limongi (2004) for high 
mass stars and van den Hoek & Groenewegen (1997) and Gavilan et al. (2005) for low 
and intermediate stars (see Fig. [3|) . The latter set of yields has now been extended to 
low metallicities ([Gavilan et al. 2 007). 

4. Results and Discussion 

The modelling result using a Scalo (1986)Scalo86 IMF is presented below in terms of 
O/H- and N/O-gradients and radial hydrogen distributions (Fig. Q] and Fig. [5]). To avoid 
over-production of oxygen the IMF was truncated at m up = 6OM for models with yields 
from van den Hoek & Groenewegen (1997) and at m up = 70M Q for models with yields 
from Gavilan et al. (2005, 2007). 

I adjusted the remaining model parameters so that the O/H-gradients were reproduced 
and a reasonable agreement with the observed hydrogen distributions were obtained. For 
the Milky Way, the properties of the solar neighbourhood were also used as constraints. 
The modelled N/O-gradients are then not in very good agreement with the observations - 
the only exceptions being NGC 598 and NGC 3031, where the models reproduce the N/O- 
gradient reasonably well, although with different sets of yields. If the abundances can be 
trusted, which is very likely, this indicates a relatively strong correlation with metallicity 
of the integrated nitrogen yields. In contemporary stellar yields for IM-stars (see Fig 1.) 
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Figure 2. Same as Fig. Q] but for NGC 2403, NGC 3031 and NGC 6946. 

there is no strong dependence on metallicity, which perhaps indicates that the IMF is 
varying. Note, also, that newer nitrogen yields, in general (see Gavilan et al. 2005)) . are 
larger than those by van den Hoek & Groenewegen (1997). 

Variations of the slope of the IMF changes the N/O-gradient, but at the same time 
the over-all oxygen abundance will change (see Fig. 2]). A Salpeter (1955) IMF, <fi(m) ~ 
mr^ x+x ' ,x — 1.35, truncated at m\ OVJ = O.O8M and m up = 5OM gives a result similar 
to that with the Scalo (1986) IMF for the O/H-gradient, but also a reasonable fit to the 
observed N/O-gradient for the Milky Way (see Fig. [4j) using the Gavilan et al. (2005, 
2007) yields. However, for the other five galaxies, this combination of IMF and yields 
does not really improve the fit to the observed abundances. Hence, resolving the O/H- 
N/O-inconsistency by changing the IMF requires that the IMF varies over the galactic 
disc, or, effectively, from galaxy to galaxy. 

A stellar-metallicity dependence in the nitrogen yields of IM stars is also a reasonable 
explanation. Since the N/O-gradients predicted by the models are too flat in three cases 
(NGC 224, NGC 2403 and NGC 6946), the yields may need to increase with metallicity. 
If the N-yields are as high as those by Gavilan et al. (2005) for IM-stars around solar 
metallicities, but as low as those by van den Hoek & Groenewegen (1997) at subsolar 
metallicities, the observed slopes might be reproduced. On the other hand, it is quite 
likely that such an ad-hoc modification of the yields may lead to gradients that are too 
steep for the Milky Way and NGC 3031. 

5. Conclusions 

Employing a standard IMF (jScalo 1986|) with a high-mass cut-off at m up — 6O/7OM 
the O/H-gradients are well-reproduced by numerical CEMs. However, the N/O-gradients 
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Figure 3. Left panels: N- and O-yields from Gavilan et al. (2005, 2007) including the previ- 
ously unpublished sets for low metallicity (left). Right panels: the same from van den Hoek & 
Groenewegen (1997). Note that the N-yields by Gavilan et al. are roughly five times larger than 
those by van den Hoek & Groenewegen. 



cannot, in general, be reproduced simultaneously. This inconsistency can be explained in 
several ways, e.g., 

(a) the production of nitrogen in intermediate-mass stars is strongly correlated with 
metallicity, 

(6) the IMF varies over the galactic disc (or, effectively, from galaxy to galaxy) such 
that the N/O-ratio is affected, 

(c) the infalling gas that forms the disc is not primordial, but significantly enriched 
with heavier elements (from Pop. Ill or Pop. II stars?), 

(d) or the P-method produces systematic errors. 

I believe that the explanations (a) and (b) are more likely then the other two. The main 
reasons are: (1) many ingredients of stellar evolution models that affects the nucleosyn- 
thetic yields (e.g., mass loss rates and the efficiency of hot-bottom burning) are still highly 
uncertain, and (2) the turn-over in the IMF at low stellar masses appears correlated with 
the Jeans mass ()van Dokkum 2008)) . (3) significant effects from chemically enriched infall 
is only expected in the outer parts of the discs, where the abundances are low, and (4) 
the P-method is in very good agreement with the standard electron-temperature method 
( jPilyugin 2001a[[Pilyugin 2001b[ ). Furthermore, the stellar yields and detailed properties 
of the IMF are probably the greatest uncertainties in CEMs in general. 
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Figure 4. Variations of the IMF. Changing the slope of a simple power-law IMF (m up = 50Mq, 
mi ow = O.O8M0 are left unchanged) will shift the O/H-gradient up or down (see left panel), 
but leave d(0/H)/dR almost unaffected. A Salpeter (1955) IMF (x — 1.35) which provides an 
acceptable fit to the observed N/O-ratios (right panel). 
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